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Abstract. In this paper we consider the classical relativistic motion of charged 
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electromagnetic knots from maps between the three-sphere and the two-sphere, we 
introduce a mean quadratic radius of the energy density distribution in order to study 
some properties of this field. We study the classical relativistic motion of electrons in 
the electromagnetic field of the Hopf map, and compute their trajectories. It is observed 
that these electrons initially at rest are strongly accelerated by the electromagnetic 
force, becoming ultrarelativistic in a period of time that depends on the knot energy 
and size. 
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1. Introduction 

As pointed out in a recent paper by Irvine and Bouwmeester p], electromagnetic knots 
are exact solutions of the classical Maxwell equations of electromagnetism in vacuum. 
They appeared firstly in a paper by Ranada [2] in 1989. Ranada himself [3] has 
used these solutions as the basic elements of a topological model of electromagnetism 
which is locally equivalent to Maxwell's standard theory but implies furthermore some 
topological quantization conditions with interesting physical meaning jH El EJ [7J E] • 

Electromagnetic knots are defined through two fundamental complex scalar fields 
(4>, 9) whose level curves coincide with the magnetic and electric lines respectively, each 
one of these lines being labelled by the constant value of the corresponding scalar. Both 
scalars are assumed to have only one value at infinity, which is equivalent to compactify 
the physical three-space to the sphere S 3 . Moreover, the complex plane is compactified 
to the sphere S 2 via stereographic projection. As a result of such compactifications, 
the scalars <fi and 9 can be interpreted, at any time, as maps S 3 — > S 2 , which can be 
classified in homotopy classes, characterized by the value of the Hopf index n [9] . It can 
be shown that the two scalars have the same Hopf index and that the magnetic and the 
electric lines are generically linked with the same Gauss linking number I. If /j, is the 
multiplicity of the level curves (i.e. the number of different magnetic or electric lines 
that have the same label or 9), then the Hopf index of both scalars is n = £fi 2 . The 
Hopf index can thus be interpreted as a generalized linking number if we define a line 
as a level curve with \i disjoint components. 

An important feature of the model is that the Faraday 2-form T = ^F^dx 11 A dx v 
and its dual *J r = | *F fMU dx fl A dx v are proportional to the two pull-backs of a, the area 
2-form in S 2 , by <p and 9, 

J 7 = —y/E(j)*a, *F = c\fa9*o, (1) 

where a is a constant introduced so that the magnetic and electric fields have correct 
dimensions and c is the velocity of light in vacuum. In the International System of 
Units, a can be expressed as a pure number times the Planck constant h~ times the 
light velocity c times the vacuum permeability [iq. As a consequence of the definitions 
(TI)), the maps and 9 are dual to one another, *(0V) = —9*a, where * is the Hodge 
or duality operator. This duality condition guarantees that both T and *J r obey the 
Maxwell equations in empty space without the need of any other requirement. 

The electromagnetic fields obtained in this way are called electromagnetic knots. 
They are radiation fields as they verify the condition E • B = 0. It can be proved 
(see [5] for the details) that any radiation field in vacuum is locally equivalent to an 
electromagnetic knot. Moreover, because of the Darboux theorem, any electromagnetic 
field in empty space can be expressed locally as the sum of two radiation fields. 
Consequently, a model of electromagnetism based on these electromagnetic knots is 
locally equivalent to Maxwell standard theory. However, its difference from the global 
point of view has interesting consequences, as are the following topological quantizations: 
(i) The electric charge of any point particle must necessarily be equal to an integer 



Motion of charged particles in an electromagnetic knot 



3 



multiple of the fundamental value g — V^ ce o ( see IS])- (ii) The electromagnetic 
helicity H = h~c(N R — N L ) is also quantized (5], where N R and N L are the classical 
expressions of the number of right- and left-handed photons contained in the field 
(i.e. N R — N L = f d 3 k(a R a R — olgil), a R (k),a L (k) being Fourier transforms of the 
vector potential in the classical theory, but creation and annihilation operator in 
the quantum version). In fact, for any electromagnetic knot, n = N R — Nl, which is a 
remarkable relation between the Hopf index (i.e. the generalized linking number) of the 
classical field and the classical case of the difference N R — Nl- (iii) The topology of the 
model implies also the quantization of the energy of the electromagnetic field in a cavity 
[?]. (iv) Electromagnetic knots are compatible with the quantization of the magnetic 
flux of a superconducting ring [8 J , which in standard theory is always an integer multiple 
of g/2, g being the value of the magnetic monopole. 

In this work we study a classical charged particle in a knotted electromagnetic 
field. We find that the particles can accelerate to light velocity. We first revise the 
construction and some physical propeties of the electromagnetic field built from the 
Hopf map between the compactified three-space (the sphere S 3 ) and the compactified 
complex plane (the sphere S 2 ). In particular, we pay attention to the electromagnetic 
energy density and how it evolves with time and introduce a mean quadratic radius 
of the energy density. Then we consider the relativistic motion of electrons in this 
electromagnetic field. We study the trajectories of the electrons and their velocities, 
showing that they become ultrarelativistic for a wide range of the electromagnetic energy 
of the knot. Finally we give some conclusions and prospects of future work. 



2. The Electromagnetic field of the Hopf fibration 

A method to find explicitly some electromagnetic knots can be found in [10J. Let 0o( r ) ; 
# (r), two complex scalar fields such that can be considered as maps (fi , 9 Q : S 3 — )■ S 2 
after identifying the physical space R 3 with S 3 and the complex plane with S 2 . They 
have to satisfy the following two conditions: 

1. The level curves of <po must be orthogonal, in each point, to the level curves of 
do, since we know that electromagnetic knots are radiation fields (E ■ B = 0). 

2. The Hopf index of (fi and of 9 are equal, H{<Pq) = H(9 ). This is necessary to 
ensure that the condition E • B = is maintained during time evolution. 

Given O an d 0o with this two conditions, we can build the magnetic and electric 
fields at t = as 

B( 0) = v 7 ^ V( ^o x V( ^o 
^ } 27rz(l + o o )2' 

EM)= 2^(iTW () 

It is convenient to work with dimensionless coordinates in the mathematical spacetime 
S 3 x R, and in S 2 . In order to do that, we define the dimensionless coordinates 
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(X,Y, Z,T), related to the physical ones (x,y,z,t) (in the SI of units that we will 
use in this work) by 

{X,Y,Z,T) = ^(x,y,z,ct), (3) 

and r 2 /Li = (x 2 + y 2 + z 2 )/Ll = X 2 + Y 2 + Z 2 = R 2 , where L is a constant with length 
dimensions. Now, let us consider the Hopf map, 

2(X + iY) 

00 = 2Z + »(iP-l) ' (4) 
whose fibres are been used basis for of knotted entanglement in reaction- 

diffusion models, in particular for a FitzHugh-Nagumo model [TTJ. We also consider the 
map corresponding to the change (X, Y, Z) y (Y, Z, X) in 

Because of their construction, it is obvious that both maps (J4]) and (j5J) have the same 
Hopf index. In fact, these maps have Hopf index n = 1 and their fibrations are mutually 
orthogonal at each point. Consequently, we can build an electromagnetic knot from these 
maps. The Cauchy data for the magnetic and electric fields are 

B(r ,0) = ( Y - XZ, -X-YZ,- l - Z * + X * + yI 



ttL 2 {1 + R*) 3 \ ' ' 2 

E(,,0) = -^^-( ^^-f-^ .-Z + XYY + XZ). (6) 



nL 2 (l + Riy 

From ([6]), two vector potentials A and C can be computed, such that B = V x A, 
E = V x C, with the results 

A(r - 0) = T M 2 if fl y ' F -- A '- 1) - 

c ( r -°) = drfW (1 '- z ' r) - (7) 

The magnetic and electric helicities of this knot are defined to be 
h m = — I A B d 3 r, 



2^o 



R? 



h e = £ -^ I C-Ed\ (8) 



2 JR3 

where Eq is the vacuum permittivity. Taking into account the Cauchy data (j6]) and the 

potentials (J7J), the electromagnetic helicity yields 

a a a . . 

h = h m + h e = — + — = —. (9) 

To find the electromagnetic knot at any time from the Cauchy data (jSJ), we use Fourier 
analysis. The fields turn out to be [TO] . 



'a 

ttLq(A 2 ~+ T 2 ) 3 



B ( r > *) = -FWTT-^ (QHi + PH 2) , 



E ^^ ^(^T 2 )3 (gH2 - PHl) ' (10) 
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where the quantities A, P, Q are defined by 

A = R2 ~^ 2 + 1 , P = T(T 2 - 3A 2 ), Q = A(A 2 - 3T 2 ), (11) 
and the vectors H4 and H 2 are 

-1-Z 2 + X 2 + (Y + T) 2 



Hi = [Y + T- XZ, -X - (Y + T)Z, 
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H 2 = ( 1 + X2 + Z ' 2 (r + T)2 - ,-Z + X(Y + T),Y + T + XZ^J . (12) 

This solution fulfills now E ■ B = and E 2 — c 2 B 2 = at any time. It is possible to 
obtain directly the electromagnetic field f lTU]) from the time-dependent expressions ([T]). 
In terms of the magnetic and the electric fields, we have 



B(r, t) = V0 ^ = / ^ w _ 96 

K,) 2vri(l + #) 2 ^ v 2mc{l + 99) 2 \dt dt 

where the time-dependent expressions of the maps and 9 are (see [5]) 
(AX - TZ) + i(AY + T(A - 1)) 
^~ (AZ + TX) + i(A(A - 1) - TF) ' 
(AF + T(A-l)) + z(AZ + TX) 

(AX — TZ) + i(A(A — 1) — TY) ' 1 ' 

Note that the level curves of both maps <fi and 9 remain linked with a Gauss linking 
number equal to 1. The evolution of the curved lines of these maps gives the evolution 
of the force lines of the magnetic and electric fields given by the expressions (TT3T) . 
The energy density of the electromagnetic field fllOp is given by, 

U(r t) - £oE2 + B2 - a ^ + X 2 + (Y + T) 2 + Z 2 ) 2 

The maximum of the energy density is located at X = Z = during time evolution. 
The function U is symmetric in the coordinates X and Z. In Figure [IJ we show some 
isosurfaces of the energy density U for times T = 0, 0.5, 1, 1.5. The energy density levels 
represented are 0.1, 0.2 and 0.3 in a/(/^ -^o) units. It can be seen how the isosurfaces 
spread as the energy density goes to zero. For time T = 1.5, the 0.3 and 0.2 levels have 
dissapeared. Note that the total electromagnetic energy of the knot 

£= f Ud 3 r = —, (16) 
J HoLo 

remains constant. At T = 0, the energy density has spherical symmetry and its 
maximum is located at the origin. As time increases, the symmetry is broken and 
the maximum is located at X = Z = and Y close to T. Approximately, the position 
of the maximum of the energy density can be found up to T = 1 at X = Z = 0, 
Y = T(l + 6T 2 )/(2 + 6T 2 ). 
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Figure 1. Evolution of three energy density levels of the electromagnetic knot of the 
Hopf fibration. From left to right and top to bottom, T = 0, T = 0.5, T = 1, T = 1.5. 
Coordinates (X, Y, Z) are dimensionless and are related to the physical coordinates 
(x, y, z) through (X,Y,Z) = (x,y,z)/L . The levels correspond to values 0.1, 0.2 and 
0.3 of the density energy in a/ipoL^) units (see Eq. (fT"5|) ). The density energy levels 
displayed are 0.1, 0.2 and 0.3 in increasing order of colour intensity. At T — 1.5, the 
levels 0.2 and 0.3 are not present. 



The energy density of the knot extends to infinity, but we can define a mean 
quadratic radius of the energy distribution as 

9 f(r — r max ) 2 U d 3 r , , 

< r >= ^jmk — • < 17) 

where r max is the position of the maximum of the distribution. At T = 0, the distribution 
has spherical symmetry and the mean quadratic radius of the distribution is given by 
r 2 > = Lq. The maximum has a value of the (dimensionless) energy density equal 
to 16/7T 2 and at a distance equal to the mean quadratic value, the dimensionless energy 
density is l/V 2 . This means that, initially, more than 70% of the energy is locallized 
inside a sphere of radius L centered at the origin. As time evolves the mean quadratic 
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radius of the distribution spreads out (its value at t = L /c or T = 1 is about 1.1 L ) 
and the position of the center is at (0,7/8,0). Note that the distribution is not well 
characterized sphere as time evolves. 

One can also compute the Poynting vector P of the electromagnetic field, obtaining 

„ f E x B o / ac \ , . 

P= / d 3 r = 0,— — ,0 . 18 

As can be seen, it has a single contribution along the y-axis. This explains why the 
maximum of the energy density moves along this axis for the electromagnetic knot 
studied in this paper. 



3. Relativistic motion of charges in the electromagnetic knot of the Hopf 
fibration 

Now we apply the electromagnetic knot studied in the previous section to the following 
situation. Suppose that this knot has been created in certain region of the space, so we 
have at t = a knot initially centered at the origin, that moves with time along the 
y-axis as we discussed in the previous section. Let us assume there are free electrons 
with negligible initial velocities in units of the light velocity c. We consider the evolution 
of these electrons under the electromagnetic knot field. 

The velocity of the electrons increases by the action of the electromagnetic field, so 
we will solve the relativistic equation for the motion of single electrons |12j , considered 
as test particles which do not affect the value of the electromagnetic field obtained from 
the Hopf fibration. The equation to be considered is 

£=-^ r ?( E +^ B -^- E >)' (19) 

where e = 1.6 x 10 -19 C is the electron charge, and m = 9.1 x 10~ 31 kg is its rest mass. 
Using dimensionless quantities and the expressions ( flOl) for the electromagnetic knot of 
the Hopf fibration, we find 

(TV _ ey/a VI - V 2 
~dT~~ nmcLo (A 2 + T 2 ) 3 

. (QH 2 - PHi + V x (QHi + PH 2 ) - V(V • (QH 2 - PH X )) . (20) 

Different possible physical situations can be studied with equation ( 120]) by changing 
the value of the dimensionless prefactor g = ey/a/ '(irmcL ). For electrons, taking into 
account that the total energy of the field is £ = 2a/(/ioA)) an d that the size of the knot 
is characterized initially by L , the prefactor in equation (|20l) can be written as 

g = ^L x0 . 15 M, (21) 
irmcLo V -^o 

where £ is measured in Joules and Lq in meters. 

In figure [2J the evolution of the position of electrons can be seen when the prefactor 

is g = 1. The initial velocity of the electrons is Vq = 0. We have taken for these 
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Figure 2. Evolution of the position of electrons under the influence of the 
electromagnetic knot of the Hopf fibration with g = 1. From left to right and top 
to bottom, T — 0.5, T — 1, T = 1.5. The energy levels shown correspond to values 
0.1, 0.2 and 0.3 of the density energy in a/^oLg) for the same times (see figure [lj. 
The initial velocity of the electrons is Vq = 0. The electron initial positions considered 
are given by R = (^,0,0), with R t = ±0.1, ±0.2, . . . , ±1, by R = (0,^,0), and by 
Ro = (0,0, Ri). After T = 1.5 these electrons moves almost freely. The final values 
of the velocities at T = 1.5 run from V m i n — 0.5300 to V max = 0.8410 in units of the 
light velocity c. 

electrons initial positions along the x-, y- and z-axis. In each axis, the dimensionless 
position is ±0.1, ±0.2, . . . ± 1. In the figure, the electron trajectories are plotted from 
T = up to T = 0.5 (in the first pannel), up to T = 1 (in the second pannel), and 
up to T = 1.5 (in the third pannel). We also include some energy density levels of the 
electromagnetic knot at the time considered, corresponding to values 0.1, 0.2 and 0.3 
in a/(/i -^o) units (see the caption of figured] for more details). By doing so, in figure [2] 
we see not only the evolution of the position of some free electrons, but also the region 
in which the interaction with the electromagnetic knot (whose center moves along the 



Motion of charged particles in an electromagnetic knot 



9 



y-axis with a velocity close to c) is more important. This fact explains some features of 
the behaviour of the electron trajectories. At the beginning, up to T = 0.5, the region 
in which the electrons are is influenced by high values of the knot energy density, so 
that the electrons are accelerated and their trajectories are curved by the influence of 
the electromagnetic force. From T = 0.5 up to T = 1, the center of the knot has moved 
along the y-axis and many electrons are now in a region in which the electromagnetic 
density is smaller. However, the electromagnetic force is still important so that many 
electrons tend to follow the motion of the center of the knot (a small fraction of the 
electrons seem to be curved towards the opposite direction). This is also the situation 
from T = 1 up to T = 1.5, during which electrons reach high values of the velocity, in 
the range from V min = 0.5300 to V max = 0.8410 in units of the light velocity c. After this 
time, the influence of the electromagnetic knot on the electron trajectories is smaller, 
and they can be considered as almost free. 

In figure [3j the evolution of the position of electrons is studied when the prefactor 
is g = 10. The initial velocity of the electrons is Vq = 0, and the initial positions in 
each axis are given by ±0.1, ±0.2, . . . ± 1 in units of the knot size L as before. As 
in figure EJ the electron trajectories are plotted in figure |3]from T = 0uptoT = 0.5 
(in the first pannel), up to T = 1 (in the second pannel), and up to T = 1.5 (in the 
third pannel), and we include some energy density levels of the electromagnetic knot at 
the time considered, corresponding to values 0.1, 0.2 and 0.3 in a/(/i -^o) units. Up to 
T = 0.5, the electrons get higher acceleration and their trajectories are more curved by 
the influence of the electromagnetic force than in the case of g — 1 (figure |2]) . From 
T = 0.5 up to T = 1, all the electrons tend to follow the motion of the center of the 
knot, a situation much clearer from T = 1 up to T = 1.5. The final velocities of these 
electrons are in the range from V min = 0.9684 to V max = 0.9942 in units of the light 
velocity c. 

In figure HI the evolution of the position of electrons can be seen when the prefactor 
is g = 100. The initial velocity and initial position of the electrons are the same as in the 
cases g = 1 (figure [2]) and g = 10 (figure [3]). As in previous cases, electron trajectories 
are plotted in figure H]from T = up to T = 0.5 (in the left-hand side pannel), up to 
T = 1 (in the center pannel), and up to T = 1.5 (in the right-hand side pannel), and we 
include some energy density levels of the electromagnetic knot at the time considered, 
corresponding to values 0.1, 0.2 and 0.3 in a/(/i £o) units. The global behaviour of the 
electron trajectories is very similar to the case in which g = 10, but now the electrons 
are even more accelerated and their trajectories are even more curved by the influence 
of the electromagnetic force. All the electrons clearly tend to follow the motion of the 
center of the knot along the y-axis, with final velocities in the range from V m i n = 0.9870 
to V m ax = 0.9999 in units of the light velocity c. 
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Figure 3. Same situation as in figure [2] with a value g = 10. From left to right and 
top to bottom, T = 0.5, T = 1, T = 1.5. The energy levels shown correspond to values 
0.1, 0.2 and 0.3 of the density energy in a/(^oLp) l° r the same times (see figure [TJ. 
The initial velocity of the electrons is Vq = 0. The electron initial positions considered 
are given by R = (^,0,0), with R t = ±0.1, ±0.2, . . . , ±1, by R = (0,^,0), and by 
Ro = (0, 0, Ri). The final values of the velocities at T = 1.5 run from V m i n — 0.9684 
to Vmax = 0.9942 in units of the light velocity c. 

4. Conclusions 

In this paper we have considered the classical relativistic motion of charged particles 
in a knotted electromagnetic field. We have seen how to construct electromagnetic 
knots from maps between the three-sphere and the two-sphere. In the particular case 
of the Hopf map, whose fibres are mutually linked, we have written the expression for 
the electromagnetic field. This is a solution of the Maxwell equations in vacuum such 
that any pair of electric lines is a link and any pair of magnetic lines is a link. We 
have considered some properties of this field, in particular the electromagnetic energy 
density. We have seen that the most part of the energy density of the knot is, at t — 0, 
localized into a sphere whose radius is the mean quadratic radius of the energy density 
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Figure 4. Same situation as in figures [2] and [3] with a value g = 100. From left to 
right and top to bottom, T = 0.5, T = 1, T = 1.5. The energy levels shown correspond 
to values 0.1, 0.2 and 0.3 of the density energy in o/(/zo£q) f° r the same times. The 
initial velocity of the electrons is Vq = 0. The electron initial positions considered are 
given by R = (Ri, 0,0), with R t = ±0.1, ±0.2, . . . , ±1, by R = (0,i? 4 ,0), and by 
Ro = (0, 0, Ri). The final values of the velocities at T = 1.5 run from V m i n — 0.9870 
to V max = 0.9999 in units of the light velocity c. 

distribution. As time evolves, the spherical symmetry is broken and the mean quadratic 
radius of the distributios spreads out. 

We have considered the relativistic motion of electrons (considered as point 
particles) in the electromagnetic field of the Hopf map. We have computed the 
trajectories of the electrons starting with zero initial velocities, and we have seen 
that these electrons are strongly accelerated by the electromagnetic force, becoming 
ultrarelativistic in a period of time that depends on the knot size. 

Finally we consider that a deeper understanding of the interaction between 
electromagnetic knots and test particles could be useful to design experiments to produce 
knots in the laboratory. 
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